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Abstract. Inhomogeneous distribution of the pinning force in superconductor re- 
sults in a magnetization asymmetry. A model considering the field distribution in 
superconductor was developed and symmetric and asymmetric magnetization 
loops of porous and textured Bii. 8 Pbo.3Sri.9Ca 2 Cu30 x were fitted. It is found that 
the thermal equilibrium magnetization realizes in crystals smaller than some size 
depending on temperature and magnetic field. 

PACS: 74.25.Ha, 74.25.Wx, 74.81.Bd 

1. Introduction 

Critical characteristics of superconductors may be derived from the dependence of magne- 
tization M on magnetic field H. However many factors affect the shape of the hysteretic M(H) 
curve. The same sample can have a various magnetization loop shape if it has been measured at 
different temperatures or up to different values of the maximal magnetic field. Also the geome- 
tric shape of the sample, its structure and homogeneity influence on the M(H) dependence [1,2]. 

At temperatures higher than ~ 10 K the high- field M(H) dependencies of high- T c super- 
conductors have a pronounced asymmetry with respect to the M = axis, so the main parts of 
hysteretic loops locate in the second and the fourth quadrants. This asymmetry grows with in- 
creasing temperature. Chen et al. [3] following approach [4] considered the total magnetization 
as the result of both the surface and the volume supercurrents, so the equilibrium magnetization 
and the surface barrier were included. As suggested in [3], the surface layer of a bulk sample, 
with a thickness of the order of the London penetration depth, is always in equilibrium. In the 
remnant volume of the sample the critical state model can be used with modified boundary con- 
ditions [3]. 
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Because the model [3] is phenomenological and has a large number of fitting parameters, 
we were interested to obtain the similar results using a different approach. Earlier, the approach 
[5] was successfully applied for symmetric hysteretic loops of T11223 and Bi2223 measured at 
low temperature [5,6]. In this article we have reexamined the model [5] and have modified it to 
fit the asymmetric magnetization loops. The resulted model explains origin of the M(H) asym- 
metry and gives results consistent with [3]. The asymmetric M(H) dependencies of the porous 
high- T c superconductor and textured one with the same chemical composition were fitted by the 
model. Our approach on the M(H) dependencies is described in Sec. 2. In Sec. 3 the magnetiza- 
tion loops of porous and textured Bii. 8 Pbo.3Sri.9Ca2Cu30 x samples are presented. Discussion and 
conclusion are given in Sec. 4. 

2. Model 

All realizing M(H) dependencies of superconductors are distributed between two extreme 
cases: a symmetric magnetization loop and a fully reversible M{H) curve. Symmetry of M(H) 
curve means that M~(H) = -Af*~(H) at high H except region near maximal applied field H m , where 
M~(H) < is the branch of magnetization curve measured during the field increasing from to 
H m , Kt(H) is the branch of magnetization curve measured during the field decreasing H m to 0. 
Symmetric M(H) dependencies are observed at low temperatures when a strong pinning of the 
vortices realizes. If there is no pinning in the sample then the M(H) dependence is fully reversi- 
ble, M~(H) = M^iH) for any H. The reversible M(H) dependencies are observed for clean super- 
conductors. High-r c superconductors demonstrate crossover from the symmetric magnetization 
curves at low 7 to the reversible M(H) dependencies at high T while pinning force decreases. 

The critical state model describes successfully symmetric M(H) dependencies in high 
magnetic fields. In version [5] it includes the next points: 

i. Magnetization of the material, which consists of a large number of cylindrical granules, 
described by the statistical ensemble (p(i?), is determined by 

oo R 

\tM(H) = -H + ( 1 - P) jU n H + 21 B% | (p{R)\ rB(r)drdR , ( 1 ) 



where P is the fraction of the material concentrated in the superconducting granules, ju n is the 
magnetic permeability of the intergranular material, cp(i?) is the distribution density of the super- 
conducting granules with respect to the radius R, B(r) is the magnetic induction in the sample. A 
homogeneous sample is described by (1) with P = and Ro = R. 

ii. The distribution of the magnetic induction B(r) in the granule is deduced from the critical 

state equation dB(r)ldr = ±(4n/c)j c (B), where j c is the critical current density. Separate sets of eq- 
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uations (2a,b,c) are resulted for three main parts of the M(H) hysteresis. 1) During the initial field 
increasing from to H m , the B(r) dependence is determined by equation: 

0(B) - 0(ji„H) = ~(4n/c) jco (R-r). (2a) 

2) For decreasing H from H m to (i.e. ht(H) branch), the B(r) dependence consists of two 
curves which are described by equations: 

0(B) - <Z>(ju n H) = (4n/c)j c0 (R-r), 

0(B) - <S>(jUnH m ) = - (4n/c)j c0 (R-r). (2b) 

3) Then, for changing H from to -H m (i.e. M(H) branch), the B(r) dependence consists of 
three parts: 

0(-B) - 0(- f i n H) = - (4n/c) jco (R-r), 
0(B) + 0(-fitf) = (4k/c) jco (R-r), 

0(B) - 0(ju n H m ) = - (4n/c)j c0 (R-r). (2c) 
Here j c o is the critical current density at H = 0, <&(B) is a nondecreasing function, such that 
d<Z>(B)/dB = (j c (B)/j c0 )- 1 [5]. 

The critical state model implies the strong pinning of vortices in a sample. In reality some 
distribution of the pinning force F p (r) is realized in granules. Therefore a region near granule 
surface can appear at high T and high H where F p is such small that vortices are not pinned. So 
the thermal equilibrium magnetization is produced by the surface supercurrent circulating in this 
region. This conclusion is coincide with idea [3,4]. To account the equilibrium magnetization we 
should modify equations 2a,b,c. We introduce new characteristic size l s determining the depth of 
surface layer where the vortices are not pinned. Boundary conditions are changed due to the sur- 
face supercurrent j s for the first equations in each set (2a), (2b) and (2c) and equation for the B(r) 
distribution in the equilibrium surface region is added: 

0(B) - 0(ju n H) = -(4n/c)j s0 (R-r), 

0(B) - 0(B S ) = -(4n/c)j c0 (R-r-l s ), (3a) 



0(B) - 0(ju n H) = -(4n/c)j s0 (R-r), 
0(B) - 0(B S ) = (4n/c)j c0 (R-r-l s ), 

0(B) - 0(jUnH m ) = - (4n/c)j c0 (R-r), (3b) 



0(-B) - 0(-n n H) = -(4tt/c) 7,0 (R-r), 
0(-B) - 0(\B S \) = - (4n/c)j c0 (R-r-ls), 
0(B) + 0(-fi„H) = (4n/c)j c0 (R-r), 

0(B) - 0(MnH m ) = - (4n/c)j c0 (R-r), (3c) 
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here j s o is the surface critical current density at H = 0, B s is the induction value at r = R-l s . The 
separate notation for j s is introduced to conformity with approach [3]. However we suppose that 
in most cases j s =j c . 

Figure 1 displays the computed B(r) dependencies after the field decreasing from H m to 
0.6H m . Curve 1 on figure 1 (solid line) is the B{r) dependence for l s = 0, curve 2 (circles) is the 
B(r) dependence for < l s < R, curve 3 on figure 1 (dotted line) is the B(r) dependence for l s > R. 
The B(r) distributions in figure 1 were calculated with the j c {B) dependence characterized by qu- 
alitatively different behaviors in scales of low and high fields [5]: 



where j c o is the critical current density at B = 0, B\ and Bo are the parameters with the induction 
measure which determine the characteristic scales. Dimensionless parameters Q > and g > 
determine the j c decrease rate. Such j c (B) dependence is convenient because it is easy integrated. 

Some M(H) dependencies computed by the obtained model are presented in figures 2 and 
3. This model is strongly adequate for the M(H) dependencies in high fields only. The surface 
barrier should be account for detailed calculation of low-field magnetization as it was realized in 
[3]. In contrast to the model [3], the presented model accounts l s , the position of the boundary 
between the equilibrium layer and the central region with strong pinning in the granules. 

The size of region, where pinning is missing, l s , depends on H and T. If the condition l s >R 
takes place for granules, the width of the magnetization loop \M~(H) - M*~(H)\ is equal to 0. The 
crossover from the symmetric magnetization to the reversible M(H) dependence can be observed 
while l s changes from to l s > R. The value of l s is believed to be about the London penetration 
depth "k [3]. Therefore if a superconducting crystal or grains in a polycrystalline sample have a 
size smaller than 2X then pinning is expected to be missing in all volume and the M(H) depen- 
dencies are reversible at any temperature. 

3. Experiment 

Synthesis of a porous Bii.gPbo.3Sri.9Ca2Cu30 x (Bi2223) is described in [7]. This material is 
formed by disordered and unpacked flakes-like crystallites with the width of ~1 |am and the leng- 
thwise sizes of -20 jam. The density of the porous Bi2223 is equal to 1 .6 g/cm . The temperature 
of superconducting transition is 1 13 K. 




(4) 
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The textured Bi2223 was prepared from the porous one [8]. The textured ceramics with the 
texture degree of 0.97 consists of ordered crystallites with the sizes of ~1 jam along the c axis 
and -10-20 jam in the ab plane. The density of textured Bi2223 is equal to 5.3 g/cm . 

Samples for magnetic measurements had cubic form with the edges of 1.6 mm. The M(H) 
dependencies were measured at temperatures of 4.2, 14, 25, 40, 60 and 80 K. The textured sam- 
ples were investigated in two orientations of field H relative to the predominant direction of the 
crystallites in the textured sample: the magnetic field is directed perpendicularly to the texturing 
plane (H\\c) and along the texturing plane (H\\ ab). 

Figure 2 demonstrates the magnetization loops of the porous Bi2223 measured at tempera- 
tures (T) of 4.2 K - 80 K. Figure 3 demonstrates the magnetization loops of the textured Bi2223 
measured at temperatures of 4.2 K - 80 K for H \\ c and H \\ ab. For the M(H) dependencies at T 
= 40, 60 and 80 K, values of the irreversibility field H in (T) are smaller than the maximal applied 
field H m = 60 kOe. At any fixed temperature, of the porous Bi2223 is smaller than the irre- 
versibility fields of the textured Bi2223 and for the textured Bi2223 H in for H || c is smaller than 
H m for H || ab. The M(H) dependencies of the textured Bi2223 are scaled according to M H ^ C (H) = 
yM H ^^ ab (yH) [9] with the anisotropy coefficient y = 2.4. Results of the detailed investigation of 
anisotropy of the texture will be presented in the other work. 

The experimental M(H) dependencies demonstrate that the asymmetry increases with tem- 
perature and field. At T= 4.2 K the magnetization loops are nearly symmetric relative to the M = 
axis at fields of 10-55 kOe, but at higher temperatures the magnetization loops are distinctly 
asymmetric. At T = 80 K the M(H) dependencies are reversible for the porous Bi2223 at fields 
higher than 1.5 kOe and for the textured Bi2223 at fields higher than 2 kOe for H || c and 4 kOe 
for H || ab. 

4. Discussion 

The experimental data were fitted in the frame of the developed model. The main fitting 
parameters are the depth of surface equilibrium layer l s and the product of j c and the average gra- 
nule radius R. We used the simple dependence l s (H) = l s {H= 0) (1 + HI(Q3H m )) that works good 
for all described curves at any T. The calculated magnetization loops depends very slightly on 
different distribution functions (p(i?) for both the textured Bi2223 and the porous Bi2223 such 
that we disregarded a dispersion of the granule parameters and used average values. 

In figures 2 and 3 the computed curves are displayed together with the experimental loops. 
There is the good agreement between the presented data at all temperatures. Some discrepancy 
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may be due to a demagnetizing factor of the sample and grains [10]. The value of the fitted l s at 
H=0 grows from /, = 200 nm at T = 4.2 K to /, ~ 3000 nm at T = 80 K (see Fig. 4a). For all 
samples, at H = 40 kOe the average radius R is smaller than values of l s at low temperatures but 
l s becomes greater than RatT above 40 K (Fig. 4b). Figure 5 shows the intragranular critical cur- 
rent density estimated from the fit with R = 10 jam for textured Bi2223 and R = 20 (im for the 
porous Bi2223. 

In conclusion, the magnetization loops of porous and textured Bi2223 were measured in 
fields up to 60 kOe at different temperatures. The M(H) dependence asymmetry increases with 
increasing temperature. To describe the asymmetric magnetization loops, the extended variant of 
the critical state model was developed. In the presented model, the asymmetry is referred to the 
single parameter l s , which determines the position of the boundary between the surface layer 
with the thermal equilibrium magnetization, and the central region with the critical state magne- 
tization. As following from the model proposed, the M(H) dependence of the sample consisting 
of sufficiently small superconducting granules should be fully reversible at all temperatures. 
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R-r [arb. units] 

Fig. 1. Dependence of induction B on depth R-r from the granule surface after field decreased 
from H m to 0.6 H m for l s = (curve 1, solid line), 0<l s <R (curve 2, circles), l s = R (curve 3, dot- 
ted line). 
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Fig. 2. Dependencies of magnetization M on magnetic field H of the porous Bi2223. Experimen- 
tal data (symbols) and computed curves (solid lines). 
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Fig. 3. Dependencies of magnetization M on magnetic field H of the textured Bi2223. Experi- 
mental data (symbols) and computed curves (solid lines). 
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Fig. 4. Dependence of boundary depth l s , divided by the average radius of granule R, on tempera- 
ture at H= (a) and H= 40 kOe (b). 
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Fig. 5. Critical current density j c versus temperature at H = (a) and H= 40 kOe (b). 
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